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We study an interacting single-level quantum dot weakly coupled to three electrodes. When two 
electrodes are biased by voltages with opposite polarities, while keeping the third lead (the stem) 
grounded, the current through the stem is a measure of electron- hole asymmetry of the dot. In this 
setup we calculate the stem current for both metallic and ferromagnetic (collinearly polarized) leads 
and discuss how the three-terminal device gives additional information compared to the usual two- 
terminal setup. We calculate both the sequential and cotunneling contribution for the currents. For 
the latter part we include a regularization procedure for the cotunneUng current, which enables us 
to also describe the behavior at the charge degeneracy points. 



I. INTRODUCTION 

A systems consisting of a quantum dot weakly cou- 
pled to external leads have been extensively studied, both 
for unpolarized^^^^^ and polarized leads. ^^^^^^^^^^^ These 
all dealt with two- terminal transport properties, with a 
third electrode acting as a capacitively connected gate. 

System with more than two terminals has a long his- 
tory in mesoscopic physics. The functionality added by 
the third lead is of importance for some applications, 
in particular the so-called Y-branch structures ^iiii^ii^ii^ 
These studies consider open mesoscopic systems, where 
correlation effects can be neglected. Some recent works 
have studied the opposite case where correlations are 
important, e.g. the study of current-current corre- 
lations in a three-terminal device, ^^^^^ Kondo peaks 
in a three/four-terminal setup ^li^ and crossed carbon 
nanotubes.^^'^^^^^ Here, we consider a multi-terminal 
quantum dot and, in particular, the non-linear response, 
which brings in new information about the structure of 
the quantum dot as well as about the magnetization of 
side branches. 

The system we study in detail consists of a single level 
quantum dot coupled to three leads. Fig. [H We allow 
collinear spin polarization of leads, so that some of them 
can be polarized while others are kept unpolarized. For 
the sake of simplicity, the polarization, if any, is assumed 
to be complete. To lift the spin degeneracy, we apply 
a magnetic field to the dot collinear with polarization 
of the leads. We study the current through the central 
junction (referred to as a stem) and its dependence on the 
applied voltage, magnetic field, position of dot's energy 
states, and orientation of leads' polarization. Except for 
the limiting case of non-interacting electrons, the exact 
solution is not known, and approximative methods are 
widely applied. In the typical experimental setup a cou- 
pling between the dot and the lead is of order of jieV 
while temperature of order of meV or higher, therefore 
it is often justified to perform perturbation expansion in 
a small parameter F/Zc^T. We cut off the perturbation 
expansion on second order. 
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FIG. 1: The three terminal device. The occupation number 
is controlled via the gate by applying the voltage Vg. The 
stem chemical potential, /Hs = 0, serves as the reference level. 
The left and right leads are biased in push-pull manner with 
/Ji — eV and /ir — —eV. 



In first order, we derive the sequential tunneling cur- 
rent through a central electrode, Eq. (|9|). When the sys- 
tem is biased symmetrically, Vi = V and Vr = —V^ the 
stem current is an even function of applied voltage and 
an odd function of the gate voltage. Fig. [H The stem 
current can be either negative or positive depending on 
the nature of the electron transport. 

To understand this, we first observe that the current 
from the left to the right lead occurs at low tempera- 
tures via only two charge states and can be predomi- 
nately either electron- or hole-like. If the energy differ- 
ence between the u -\- 1 and u electron states is positive, 
the transport is electron-like and hole-like if it is nega- 
tive, with u standing for the occupation number. For the 
electron-like case the state can decay via an electron 
leaving through the stem, whereas for the hole-like case, 
the ly state can decay by the electron entering from the 
stem. Hence the stem current is zero for the electron-hole 
symmetric case. 

The three terminal setup therefore measures the 
electron-hole asymmetry, and in this respect it is similar 
to thermopower.^^ Interestingly, a change between the 
two situations above, can be induced by the magnetic 
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FIG. 2: (Color online) The stem current (solid, red) being the 
sum of the sequential part (dotted, blue) and the cotunneling 
(dashed, green). The sequential contribution gives sufficient 
qualitative description. For this plot leads are unpolarized 
with Ti ^Tr = Ts = r/3, U 20r, B 5r, ksT 2r and 

v/r = 1. 

field applied to the dot which modifies the physics signif- 
icantly, see Fig. [5l Moreover, in a three-terminal setup 
with magnetized side branches and the non-magnetic 
stem the current value enables one to distinguish among 
four magnetization alignments. Fig. [6l 

In second order, the total current is increased by vir- 
tual cotunneling processes for these values of the gate 
voltage for which the sequential tunneling is exponen- 
tially suppressed. Furthermore, the cotunneling correc- 
tion lowers the sequential current maxima and gives rise 
to additional broadening of order of the sum of all cou- 
plings Fg +r/ + r^, and hence the total broadening is the 
sum of the thermal part and the latter. However, we will 
show that the cotunneling current does not change the 
sensitivity towards electron-hole asymmetry, and there- 
fore the qualitative picture described above still holds to 
higher order in tunneling, see Fig. O 

In the Appendix the exact solution for non-interacting 
electrons is compared with the perturbative treatment 
and even though the renormalization of level positions 
also contributing to the second order current was ne- 
glected, we achieve good agreement between the two. 

II. MODEL HAMILTONIAN 

The model Hamiltonian of the quantum dot connected 
to any number of unpolarized leads is 

H = Hl^Ht^Hd, (1) 

where 

= E = E E ^JMcU„Cj,ka (2) 

j j ka 

describes the uncoupled leads indexed by j while 
and Cj^ka form an orthogonal set of creation and annihi- 



lation operators in lead j. 

The dot is described by the Hamiltonian 

Hd = ^(eo - (jB)dld^ + Un^ni (3) 

with a single electron of the orbital energy cq, the 
Coulomb repulsion energy [/, magnetic field a = 1 
for a spin- up electron and cr = —1 for a spin-down elec- 
tron, (ij. and da- forming a set of orthogonal creation and 
annihilation operators for the dot and Ua- = d^^d^ being 
an occupation number operator. The tunneling processes 
between the leads and the central region are taken into 
account via the tunneling Hamiltonian 

= E J2(hMc],kad. + h.c), (4) 

j ka 

where tj^ka is a spin-dependent tunneling amplitude. 

To include the complete magnetization of one of the 
leads, say lead j, one should limit the corresponding 
Hilbert space to one spin direction replacing a by | or |. 



III. SEQUENTIAL TUNNELING REGIME 

The simplest situation we study is the sequential tun- 
neling regime, also called the weak tunneling regime. It 
is assumed that the time between tunneling events is the 
largest time scale in the problem, so that there is no co- 
herence between successive tunneling processes. If there 
is no bias applied, the distribution function of different 
states is given by the equilibrium Gibbs function. With 
an applied voltage difference between the electrodes the 
induced non-equilibrium distribution function needs to 
be determined. 

To this end, we calculate the transition rates between 
different dot's states. Since we consider the weak tun- 
neling regime, the Fermi's Golden Rule is sufficient to 
tackle this problem. We define the transition rates F^^ 
as the rate for a process that changes the state of the dot 
from V io ji due to tunneling through junction j. Assum- 
ing a continuous density of states inside the electrodes it 
follows that the tunneling rates are proportionate to the 
Fermi function, F^+i^j, = r^n(ej,+i — — fij), if an elec- 
tron jumps onto the dot, and r^_i^^^ = r^(l — n{eiy-i — 

— /ij)) for the opposite process. Here fij is a chemical 
potential of reservoir j, F^ = 27r|tjppj with pi being the 
density of states in lead j. To include the collinear polar- 
izations Pi of the leads spin-dependent tunneling rates 
^ia = ^r?(l + aPi) are introduced. 

Having found the transitions rates, we can write down 
master equations describing the dynamical behavior of 
the distribution function Pj^ and since we are only inter- 
ested in the steady state solution, we have 
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where F^j, = F^^ is a total contribution from all 
the leads to the transition rate from state u to state /i, 
while Po, P|, P|, P2 are the empty, spin- up, spin-down 
and double occupation state distribution functions re- 
spectively. The terms with a minus sign give the rate 
at which a given state on the left-hand side decays while 
the plus-sign terms describe the opposite processes. In 
addition to these equations, the probability conservation 
law Po + P| + P| + P2 = 1 has to be used. 

The knowledge of the distribution functions allows one 
to find the net current fiowing through junction i 

^i = -|E(rt+i,.-rt-i,JP.- (6) 

u 

This should be understood as a difference between the 
number of electrons incoming to the dot and deoccupat- 
ing it, times the electron charge, — e < 0. 

We choose the Fermi level of the stem reservoir to be 
our reference level, jUg = 0, throughout the text and in all 
diagrams. 

The solutions to Eqs. (O are 

Po = -^(FotFt2(Fot + F2J + Fo|F|2(FoT +r2T)), (7a) 

n = ;^(rToroT(r|2 + f^2) + r^2F2T(F|o + r^o)), (rb) 

Pi = ^ (F^oFoT (F|2 + Fr2 ) + F12F2T (F|o + r^o)) , (7c) 
P2 = -^(F2TrTo(FoT + F2J + F2|F|o(FoT +r2T)), (7d) 
with 

^ = (r^oFoT + F2tF|o)(F|2 + F21) 
+ (roTF|2 + F|2F2t)(F|o + Fot) 
+ (FotF|2 + F|2F2t)(F|o + Fo|) 
+ (F|oFoT + F2|F|o)(F|2 + F2t). (8) 

Using Eq. (|6|) the stem current becomes 

+{Tli-nOPi-{ri2 + ^2)P2)- (9) 

The above expression is a complicated combination of 
the Fermi functions and analytic treatment is not further 
possible. We discuss numerical results in Sec. |Vl 



IV. COTUNNELING REGIME 

As discussed in the introduction, one could expect 
the cotunneling current to modify the lowest order cur- 
rent significantly. However, this turns out not to be the 
case, since the cotunneling current also has nodes in the 
particle-hole symmetric points and only give small cor- 
rections to the overall shape of the current versus gate or 
voltage curves. 

In the cotunneling regime two electron processes come 
into play. An electron is transferred between lead j and 
lead i (via an intermediate classically forbidden state) in 
two successive tunneling events, across the quantum dot. 
For the calculation of the two-electron rates, we need 
to consider the different starting configurations (empty, 
single- and double-occupied). The probabilities that 
a given state is occupied are given by Eqs. ([7|). Fur- 
thermore, there are two types of processes that should 
be considered, i.e. non-spin-fiip processes that do not 
change dot's magnetization, and spin-fiip processes lead- 
ing to reversal of spin direction, see Fig. [3l The latter 
contribute directly to the current as well as modify the 
probabilities via spin-fiips caused by the interaction 
of the dot with the lead. We find the rates for these pro- 
cesses employing the generalized Fermi's Golden Rule 

rfi = -\{f\T\i)\'s{Ef-E,) (10) 

for the transition from the initial state \i) of energy Ei to 
the final state |/) of energy Ef. The transition operator 
T is defined as 

T = Ht + HtjAjtT (11) 

with Hq = Hl + Hd- The first non- vanishing transition 
rate for the process which transfers the electron between 
the electrodes appears in the perturbation expansion in 
second order, thus the cotunneling process is quadratic 
in couplings, F. 

The cotunneling events are two-particle processes and 
therefore occur only between pairs of the leads. Hence, 
the cotunneling current through junction i can be ex- 
pressed as the sum of the currents between any pair 
of leads J^^^ i.e. the currents between lead i and j in 
regime z/, weighted by the appropriate probability Pjj 

Ji = E(J2P'^^^ + T.P''W')- (12) 

j cr 

The sum runs over dot's states u = 0,t,i,2, and for 
emphasis we divided the current into the non-spin-fiip 
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part J^ - and the spin-flip part J^-^^ (present for the single 
occupation only). 



A. Non-spin-flip cotunneling current 

We exemplify the derivation of the current between 
a pair of leads by the case of the empty dot. The initial 
state \i) = . . . , i^AT, 0) consists of a tensor product of 
lead's states and the dot's state |0). The electron can 
be transferred from lead i into lead j via either a spin-up 
or spin-down state of the dot, depending on a spin of the 
electron entering the intermediate region. Because the 
two corresponding flnal states of the leads are different 

I/) = |z^l,...,Z^i-Cr, ...,Zy^-+Cr, ...,i/7V,0) = Ci,/cac]^fc,^|i), 

there is no interference between electron's paths and we 
flnd the rates F^^^^ for these two processes separately 
and add them up to get the total tunneling rate in this 
regime. Thus, the current between the pair of leads i, j 
is given by 



-fij — 



(13) 



that is the rate for the process bringing the electron 
from lead i to lead j through the empty state, minus the 
rate for the opposite process, multiplied by the electron 
charge, — e. We substitute the tunneling Hamiltonian 
Ht^ Eq. (HI), into {f\T\i) and after a number of stan- 
dard calculations, where we use (i) that the distribution 
function of electron states in the leads are given by the 
Fermi-Dirac distribution functions and (ii) the assump- 
tion of a constant tunneling density of states, we arrive 
at 



1. 



-pu — F^ 



1 



(14) 



Employing Eq. (p!3|) . the current between lead i and j 
through the empty dot becomes 



70. - -::r' 



J — ( 
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1 



X (^(^ - Mi) - ^(^ - Mj)), (15a) 



where we took F^ = F^^ = F^^ . If one neglects the inelas- 
tic processes, the currents through the single occupied 
dot may be found by analogy 



Tcr _ __-rO-pO 



d^ 



with a = — cr. The expression for the current in the 
presence of spin-flipping will be derived below as it needs 
more attention. 
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FIG. 3: The comparison of a cotunneling event for the dot 
occupied by a spin-up electron with the spin-flip process (on 
the right), and without (on the left). In case of inelastic 
processes, the final state of the dot and the leads is the same 
for both paths resulting in the interference. 



Finally, the currents Jf, in the remaining regime of the 



double occupation reads 



t2 _ _0-pO-pO 



di 



These expressions are, however, divergent and we need 
to improve on the second order perturbation theory (in 
F) to regularize these divergences. Before we follow the 
regularization procedure in Sec. \IV C\ we study the cur- 
rent flowing with reversal of the spin accumulated on the 
dot. 



B. Spin-flip cotunneling current 

As long as we consider the limit of non-interacting elec- 
trons, = 0, the inelastic processes do not affect the 
current at all, and the above description is completely 
sufficient. In the Appendix we compare the perturbative 
result with the exact results obtainable when neglecting 
interactions. 

The derivation of the tunneling rate and the current 
between a pair of leads in the presence of spin-ffip pro- 
cesses is similar to one used in the previous section. Here, 
we just cite the result 



j^^f = -e(r 



_ r 



(16) 



with the tunnehng rate 



1 r° 
h J_, 



1 



(17) 



This expression differs from that for the cotunneling 
current without spin-flip processes, Eq. (|15bp . as the 
interference term appears there. Fig. [3l Furthermore, 
the chemical potential of lead j is shifted due to the 
energy gap, A^^r = — resulting from the Zee- 
man splitting. It is convenient to rewrite the expression 
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(n — /i j — A^^cr ) — ^ — Mi ) ) so that it has the same form 
as Eqs. ([15]), ie. it is proportional to the Fermi functions' 
difference, with ub standing for the Boltzmann function. 
In the next section it wih come to Ught that terms of this 
type undergo the regularization scheme. 



C. Regularization procedure 

In general for the non-spin- flip cotunneling current, the 
problematic integrand is a product of the divergent term 
(C ~ Co)~^ and the Fermi functions' difference, denoted 
as f(i) = n{i - iii) - n{i - iij). 

To deal with this divergence we follow the regulariza- 
tion scheme proposed by Turek and Matveev^^^^^ and add 
to the denominator a life-time broadening, rf ^ describ- 
ing the tunneling broadening of the intermediate state, 
i.e. T] oc F. The whole trick is to partition this inte- 
gral into two parts, from which the first can be a pos- 
teriori identified as the energy conserving process, and 
hence, the sequential tunneling contribution while the 
second term describes the regularized cotunneling pro- 
cesses. This is done as 



/ 



(?-Co)2 



/(Co) 



- Co) 

r/(6) 



no - mo) 



lim 



fio - mo) 

(e-?o)2 + r?2 



(18) 



where the last line is in the limit of small rj. The first 
term corresponds to the transitions on resonance (or 
"on-shell"), where energy is conserved in each tunnel- 
ing event. These processes give rise to a current which 
is linear in F after inserting into Eq. ([12]) and using that 
77 oc F. We can therefore omit the first term, since it has 
already been calculated within the much simpler master 
equation scheme in Sec. I/I/I The remaining term now 
corresponds to the proper regularized second order con- 
tribution to the current. At this point it is important to 
realize that this procedure does not capture all second 
order terms, because the renormalization of the dot state 
due to tunneling is not included^i The effect of renormal- 
ization, however, is easily incorporated by adding these 
terms (linear in F) to the energies of the dot states when 
doing the master equations. The second order correction 
(in F) can then be extracted from the master equation 
result. Here we do not incorporated them since they are 
unimportant and merely give a shift of the gate voltage. 
On the other hand, for some cases, e.g. non-collinear 
magnetization, they cannot be neglected because they 
give rise to off-diagonal elements in the density matrix.^ 
Now returning to the expression ([T8|) we can evaluate 



it analytically by employing the useful identity 



v^+J "^{^-^0)^ + ^' V^O+d^o 



(19) 

and summing over the residues. The final result becomes 



.(e-eo)2 



Re 



^(*i(eo,Mi)-*i(?o, 



Mi, 



(20) 



The arrow indicates that the divergent integral has been 
regularized by the procedure explained above. For con- 
venience we introduced a shorthand notation for the 
polygamma function of n-th order 

(21) 

Below we list the non-spin-flip currents, Eqs. ([T5|) after 
the regularization (for more detailed derivation of regu- 
larized formulas see Ref. 24) 

+ ^(*i(ei,rt)-*i(ei,Aii))], (22a) 

J!^. = - £r°r°Re[A(viy,(e,,^,) _ M/i(e.,^,.)) 

+ + U,lii) - 4-1 (e. + U,^ij))], (22b) 

J2 = - lrOrORe[^(M/i(eT + U,^,) - M/i(eT + U,^,)) 

+ ^ (*i (q + U, ^li) - vE-i (e^ + [/, Ai,))] . (22c) 

The current in the presence of the inelastic scatter- 
ing, Eq. ([T6|) . causes more problems as the interference 
term emerges. Using the partial fraction decomposition 
it turns to be 

2 
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1 



^-6 ^-6 
1 



(e-6 



(23) 

with the new type of divergence 

(^-6)"^. The identity 

ii-^l)^+V^ ri^O+ J 



lim / 



IT] 

where ^(^) = (^ — 6 + j~z^)/(05 allows one to find the 
regularized expression for this type of divergence as well 
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o((6,/i,)-^o(6,/^z))] (25) 
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FIG. 4: The stem current in function of the bias voltage V 
and the gate voltage eo. The leads are unpolarized and equally 
coupled to the dot, with coupUng's strength r/3. The picture 
without the magnetic field (on the left) differs from this with 
field B — 5r (on the right). The interaction energy U — 20r 
and temperature ksT = F. The grey scale in the far right 
describes the stem current value in units eV/h. 



and after some algebra the inelastic tunneling rates be- 
come 



1 



xRe 



+ A^^) - ^i(e^,/ii)) 

2 

-- (^o(ea + U, iij + A^^) - ^o(ea + U, in)) 



--(^o(ea,/ii + A^^) -^o(ea,/ii)) 



(26) 



which substituted to Eq. (p!6|) gives a proper limit of non- 
interacting electrons, limt/^o Jij^^ = 0- Note, that cur- 
rent between each set of leads (involving spin's reversal 
or not) fulfils the principle of detailed balance (vanishes 
for jiii = jij) and, hence, the total cotunneling current is 
equal to zero, when the chemical potentials of the leads 
are at the same level. 



V. RESULTS 

In Fig. [4] the stem current dependence upon the bias 
voltage V and the gate voltage eo for no magnetic field 
(the left plot) and for magnetic field B = bV (right) is 
shown. Bright regions correspond to positive values of 
the current while dark areas to negative ones. It is ap- 
parent that the stem current is an even function of bias V 
and an odd function of the orbital energy eo with respect 
to the particle- hole symmetry line eo = —U/2. Therefore, 
we restrict our discussion to the parts of the plots where 
y > and eo > — 1//2. By performing the particle- hole 
transformation, the behavior in region where eo > —U/2 
holds, can be mapped onto the remaining area of the 
y — eo space. 

Consider the case without the magnetic field first. If 
eo < 0, the dot is in the single occupied state, and un- 
til V < the transport is blocked because electrons 
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FIG. 5: (Color online) The stem current for different values 
of magnetic field B is plotted. Interestingly, there is a sign 
change in the vicinity of 1/ = for B — cq. See text for 
explanation. The gate voltage is eo = 6F. Other parameters 
are as in the previous figure. 



from the left branch do not have enough energy to over- 
come the Coulomb blockade. Increasing voltage makes 
the transport out of the stem possible, since electrons 
may escape to the empty states in the right lead. How- 
ever, further increase of the voltage stops the stem cur- 
rent again, because for V > eo + (when two electrons 
excitations become possible) the left reservoir supplies 
the stem with electrons while at the same time the elec- 
trons from the stem move into the right reservoir. These 
currents cancel one another and there is no net stem cur- 
rent. The situation differs for eo > 0. The current in- 
creases once V > eo, when dot is excited to the single 
occupied state, and again when F > eo + t/, where the 
electrons transverse through the double occupied state. 
In both cases the direction of the flux is directed into the 
stem and carried by the left-lead electrons. Clearly, for 
the negative voltages the role of the right and left reser- 
voir interchanges, but nevertheless the direction of the 
stem current is unaffected. 

In the presence of the magnetic field more complex 
structure emerges due to two new excitations coming 
into play (as the single occupied state is no longer de- 
generate). It is convenient to divide the right part of the 
diagram into four vertical strips eo G [—U/2^—U/2-\-B]^ 
[-U/2 + 5, 0], [0, B] and eo > B. 

In the first region, the negative stem current flows in 
the narrow range of voltages. Below the lower threshold 

V < — eo + 5, the spin- up electron occupying the dot 
stops the current until the required state in the right 
reservoir becomes available. Then, both spin-up and 
spin-down electrons can participate in transport. For 

V > €o -\- U -\- B the electrons from the left branch enter 
and compensate the current carried from the stem to the 
right lead and, hence, the net stem current vanishes. 

In the second region, from V = OtoV = — eo + 5 
the dot is occupied by the spin-up electron which blocks 
the current. Beyond V = —eo -\- B the spin- up states in 
the right lead become available and this electron moves 
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FIG. 6: The stem current for different polarizations of leads 
and magnetic field 5 = SF. The central branch is always 
unpolarized. In the upper row, both left and right leads are 
spin-up polarized, Pi = Pr = 1 (the left plot), or spin-down 
polarized. Pi = Pr = —1 (the right plot). Due to nonzero 
magnetic field, the spin symmetry is broken and the stem 
current is no longer an odd function of the gate voltage with 
respect to eo = —U/2. In the lower row, polarization of the 
left and right leads is antiparallel {Pi = 1, Pr = — 1 on the left, 
and = —1, Pr = 1 on the right). In this case, the current 
symmetry in V is violated, because the dot couples to the left 
and right leads asymmetrically. The other parameters are as 
in the previous figure. 



towards the right reservoir. Further increase of the bias 
voltage results in a small increase of the current at V = 
eo^U — B, when spin- up electrons from the left electrode 
can pass through the Coulomb blockade. Finally, there 
is a sign change of the stem current at F > eo + /7 + 5, 
because electrons of both spin directions go into the stem 
from the left lead while still only spin-up electrons can 
move into the right reservoir. 

The negative stem current in the third strip starts to 
flow when V > —eo-\-B due to the spin- up electrons 
moving out of the stem into the right reservoir. This 
current is risen by the spin-down electrons from the left 
lead at V = eo + 5 and also by electrons moving through 
the double occupied state (for both V = cq^U — B and 

V = eo^U^B). 

Eventually, in the last strip the stem current has four 
steps. These are for V = cq — B, V = cq -\- B when the 
excitations of the empty dot to spin-up and spin-down 
state respectively become energetically allowed, and for 

V = eo^U — B,V = eo-\-U-\-B when electrons from the 
right lead have enough energy to overcome the Coulomb 
blockade (due to the spin-down and then spin-up electron 
on the dot). 

Fig. [5] shows how the current is affected by the mag- 
netic field applied to the dot. The stem current remains 
positive as long as 5 < eo for any bias. The character 



of the function changes when B crosses eo- This corre- 
sponds to the situation when the single occupancy state 
has the same energy as the empty state (the electron-hole 
symmetry point). Tuning the magnetic field to 5 = eo 
enables one to measure directly electron-hole transport 
fluctuations. 

Providing B > cq^ the spin- up state of the dot is below 
the Fermi level of the stem and for small voltages the 
negative current is most likely to occur. However, the in- 
creasing voltage forces more and more electrons to move 
out of the left lead into the central branch. These elec- 
trons overbalance those heading in the opposite direction 
and the current sign changes - it is apparent for B = 8T. 

In two-lead systems the tunneling magnetoresistance 
(TMR) have widely been studied. ^^^^^ Defining a sim- 
ilar quantity in the multiprobe setup is rather vague. 
Nonetheless, it turns out that the current through the 
central, unpolarized lead has interesting features while 
varying magnetization of the other electrodes. The re- 
sults are shown in Fig. [6l 

The parallel magnetization (the upper part) destroys 
the antisymmetry in the gate voltage whereas the an- 
tiparallel alignment (the lower part) breaks the bias volt- 
age symmetry. It is worth to note that the third lead 
gives not only information about the relative magnetiza- 
tion of two remaining electrodes, but also allows one to 
determine polarization of each lead. This means that in 
that kind of a device we can switch among four different 
states and directly read out the information decoded in 
the two-lead system with the third probe. Experimen- 
tally it might be more convenient to swap the magnetic 
field on the dot rather than to change polarization of the 
side leads. 



VI. CONCLUSIONS 

In conclusion, we have studied Coulomb blockade in 
a three terminal device. We have focused on the bias sit- 
uation where two leads have opposite voltage, while the 
third lead (the stem) is grounded. This setup allows for 
a direct measure of the electron-hole asymmetry of the 
quantum dot system and we have made detailed calcula- 
tions of the stem current in the gate voltage-bias voltage 
plane. We have pointed to a number of predictions, that 
can be experimentally tested. 

Furthermore, we have considered the spin polarized 
case, where a number of detailed experimental proposals 
have been presented. In particular, a large difference 
in predicted pattern between parallel and anti-parallel 
configurations is seen. 

Finally, we have checked, using a regularized version 
of the usual cotunneling formalism, that our predictions 
are not significantly altered by cotunneling corrections. 
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FIG. 7: (Color online) The strict result for the stem current 
(thick, black) compared to the perturbative result: sequential 
(dotted, blue), cotunneling (dashed, green), and their sum 
(thin, red). Parameters are as in Fig. (|2]) (except for [7 = 0). 



APPENDIX: NON-INTERACTING ELECTRONS 
LIMIT 

The limiting case of non-interacting electrons gives the 
opportunity to compare the results we get within the 
framework of the second order perturbation theory with 
the exact result containing higher order terms in cou- 
plings. The quality of the cut-off on the second order 
terms can be estimated. We derive the generalization of 
the Meir-Wingreen formula that will include the indefi- 
nite number of leads and calculate the current Ji through 
one of the junctions. 

The main result of Meir and Wingreen^^ is the current 
through lead i in the presence of interactions 



J-- 

(A.l) 

This can be generalized to the many lead systems em- 
ploying the obvious identity = jJi — (1 — 7) Xl^^ 
following from Kirchhoff's law. Providing all the cou- 
plings are proportionate, that is = XijTj where Xij 
are constants, one may eliminate the lesser Green func- 
tion G< by proper selection of 7 showing that 



Ji 



rr 



(A.2) 

where we introduced a shorthand notation F = F^ and 
wrote the formula for the current Ji in terms of a spectral 



function A = i(G^-G^). 
For non-interacting electrons 



(A.3) 



describes a spin dependent spectral function, which sub- 
stituted to the generalized Meir-Wingreen formula gives 
the exact current in lead i for the multiterminal device, 
valid to any order in couplings 



Ji 



h 



(A.4) 

We note in passing that the same result comes from the 
non-interacting Landauer-Biitticker formalism for any 
number of electrodes. 

Using the similar methods as described in Sec. IIVC[ 
we derive the current through junction i that depends on 
coupling F up to any order 



and 



*o(e<T, Mi, r) - *o(e<T, Mi> T) 
(A.5) 



(A.6) 



In Fig. [71 the exact result is compared with the per- 
turbative one. The good agreement between these two 
justifies the choice of the method, and reassures that the 
second order perturbation theory gives not only excel- 
lent quantitative description, reachable also within the 
sequential tunneling framework, but also good qualitative 
estimation, at least for non-interacting electrons. Note, 
however, that this conclusion is only valid in the regime 
where F <C ksT^ whereas in the opposite limit the per- 
turbation expansion clearly fails, and one does not expect 
the good agreement. 
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